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Abstract 

When  operating  the  polymer  electrolyte  membrane  (PEM)  fuel  cell  stack,  temperatures  in  the  stack  continuously  change  as  the  load  current 
varies.  The  temperature  directly  affects  the  rate  of  chemical  reactions  and  transport  of  water  and  reactants.  Elevated  temperature  increases  the 
mobility  of  water  vapor,  which  reduces  the  ohmic  over-potential  in  the  membrane  and  eases  removal  of  water  produced.  Adversely,  the  high 
temperature  might  impose  thermal  stress  on  the  membrane  and  cathode  catalyst  and  cause  degradation.  Conversely,  excessive  supply  of  coolants 
lowers  the  temperature  in  the  stack  and  reduces  the  rate  of  the  chemical  reactions  and  water  activity.  Corresponding  parasitic  power  dissipated  at 
the  electrical  coolant  pump  increases  and  overall  efficiency  of  the  power  system  drops. 

Therefore,  proper  design  of  a  control  for  the  coolant  flow  plays  an  important  role  in  ensuring  highly  reliable  and  efficient  operations  of  the  fuel 
cell  system. 

Herein,  we  propose  a  new  temperature  control  strategy  based  on  a  thermal  circuit.  The  proposed  thermal  circuit  consists  of  a  bypass  valve, 
a  radiator  with  a  fan,  a  reservoir  and  a  coolant  pump,  while  a  blower  and  inlet  and  outlet  manifolds  are  components  of  the  air  supply  system. 
Classic  proportional  and  integral  (PI)  controllers  and  a  state  feedback  control  for  the  thermal  circuit  were  used  in  the  design.  In  addition,  the  heat 
source  term,  which  is  dependent  upon  the  load  current,  was  feed-forwarded  to  the  closed  loop  and  the  temperature  effects  on  the  air  flow  rate  were 
minimized. 

The  dynamics  and  performance  of  the  designed  controllers  were  evaluated  and  analyzed  by  computer  simulations  using  developed  dynamic 
fuel  cell  system  models,  where  a  multi-step  current  and  an  experimental  current  profile  measured  at  the  federal  urban  driving  schedule  (FUDS) 
were  applied.  The  results  show  that  the  proposed  control  strategy  cannot  only  suppress  a  temperature  rise  in  the  catalyst  layer  and  prevent  oxygen 
starvation,  but  also  reduce  the  parasitic  power  dissipated  for  operation  of  the  air  blower  and  coolant  pump. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  PEM  fuel  cell  is  a  device  used  to  generate  electrical 
power  via  a  chemical  reaction  whose  only  byproducts  are  heat 
and  water.  Thus,  the  fuel  cell  does  not  produce  air  pollution,  an 
advantage  that  might  allow  the  PEM  fuel  cell  to  be  considered  as 
a  potential  alternative  energy  source  for  future  automotive  and 
stationary  applications. 

To  replace  conventional  power  sources,  fuel  cell  systems 
must  exhibit  exceptional  performance,  efficiency  and  reliabil¬ 
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ity.  Important  properties  that  affect  the  performance  of  PEM 
fuel  cells  include  transport  of  reactants  and  byproducts,  rejec¬ 
tion  of  the  heat  generated  by  electrochemical  reactions  or  by  the 
current  passing  through  the  fuel  cell  and  control  of  humidity  to 
maintain  adequate  electrolyte  conductivity  [1]. 

Balance-of-plant  (BOP)  is  a  group  of  system  components 
that  supply  reactants,  remove  the  heat  generated,  remove  the 
water  produced  and  control  the  actuators.  Typical  components 
for  operation  of  a  PEM  fuel  cell  system  include  a  hydrogen 
tank  to  store  fuel,  an  air  compressor  or  blower  along  with  the 
inlet/outlet  manifolds,  a  humidifier  to  supply  humidified  air,  a 
bypass  valve,  a  radiator  with  a  fan,  a  reservoir,  a  coolant  pump, 
several  control  valves  and  controllers  to  reject  the  heat.  A  typical 
configuration  is  shown  in  Fig.  1. 

As  matter  of  fact,  the  temperature  effects  on  performance 
and  reliability  are  one  of  important  issues.  When  chemical 
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Nomenclature 

A  area  (m2) 

b  membrane  extension  coefficient 

C  mas  s  concentration  (kg  m-  3 ) 

Cp  specific  heat  ( J  kg- 1  K- 1 ) 

F  Faraday  number  (A  s  mol- 1 ) 

Fr  radiator  frontal  area  (m2) 

i  current  density  (A  m-2) 

h  heat  transfer  coefficient  (W  m-2  K- 1 ) 

J  rotational  inertia  (kg  m2) 

m  mass  (kg) 

M  molar  mass  (kg  mol  - 1 ) 

n  number 

N  mole  flux  (mol  s-1  m-2) 

p  (partial)  pressure  (Pa) 

Q  heat  transfer  (W) 

R  Universal  gas  constant  (J  kg- 1  K- 1 ) 

R  resistance  (Q) 

s  entropy  ( J  mol- 1  K- 1 ) 

t  thickness  (m) 

T  temperature  (K) 

W  mass  flow  (kg  s-1) 

Subscripts 
amb  ambient 

an  anode 

bl  blower 

c  coolant 

ca  cathode 

catl  catalyst  layer 

cond  conduction 

conv  convection 

cv  control  volume 

diff  diffusion 

ele  electro 

g  gas 

i  index 

m  motor 

membr  membrane  layer 
rad  radiator 

res  reservoir 

s  stator 

sou  source 

st  stack 

Greek  symbols 
s  porosity 

y)  efficiency 

X  water  content  (ratio) 

p  density  (kg  m- 3 ) 

r  tortuosity,  torque  (N  m) 

ip  flux  linkage  (V  s  rad- 1 ) 

co  angular  velocity  (rad  s - 1 ) 


reactions  occur  and  transport  of  charges  begins,  heat  is  gener¬ 
ated  in  the  cells  and  temperature  rises.  The  elevated  temperature 
increases  the  rate  of  chemical  reaction  and  water  activity  that 
affects  transport  of  charges  and  reactants.  However,  the  tem¬ 
perature  should  be  limited  within  a  working  temperature  that 
does  not  deteriorate  material  properties  of  components.  When 
loads  for  the  stack  continuously  change,  the  temperature  pro¬ 
file  in  the  stack  varies  and  the  limit  set  can  be  exceeded.  The 
importance  of  the  temperature  dynamics  was  recognized  and 
studied  by  some  authors  [2,3]  in  addition  to  the  temperature 
effects  on  reliability  of  membranes  and  catalysts.  A  high  work¬ 
ing  temperature  at  a  given  current  load  increases  water  flux  in 
the  membrane  from  the  anode  to  the  cathode  and  likely  causes 
dehydration  [4].  The  rate  of  Pt  particle  size  growth  gets  larger 
at  elevated  temperature  and  resulting  loss  of  platinum  surface 
area  causes  efficiency  degradation  [5].  Therefore,  management 
of  temperature  by  controlling  the  coolant  flow  plays  an  impor¬ 
tant  role  in  ensuring  highly  reliable  and  efficient  operations 
of  the  fuel  cell  system.  We  have  attempted  to  find  a  controls 
law  that  allows  for  securing  safe  and  efficient  operation  of  a 
stack. 

2.  Models  for  a  PEM  fuel  cell  system 

Control  of  a  fuel  cell  power  system  requires  a  better 
understanding  of  the  dynamic  behavior  of  the  stack  that 
interacts  with  different  BOP  components.  Due  to  the  com¬ 
plexity  of  the  system,  dynamic  models  are  used  to  efficiently 
design  and  effectively  assess  controllers.  Models  describing 
the  dynamic  behavior  of  the  PEM  fuel  cell  stack,  air  supply 
and  thermal  systems  are  briefly  described  in  the  following  sec¬ 
tions. 

2.7.  Fuel  cell  stack 

Most  fuel  cell  models  describing  physical  behaviors  of  a 
PEM  fuel  cell  are  based  on  either  empirical  equations  whose 
parameters  are  obtained  by  fitting  them  to  the  curve  of  a  spe¬ 
cific  polarization  characteristic  [6]  or  the  computational  fluid 
dynamics  (CFD)  that  basically  solves  mass  and  charge  trans¬ 
port  in  a  cell  described  by  using  the  Navier-Stokes  equations 
[7-9].  The  former  has  been  used  to  design  a  controller  for  the 
air  supply  system  of  PEM  fuel  cell.  The  lacking  dynamics  of 
a  fuel  cell  can  be  improved  by  reflecting  the  charging  and  dis¬ 
charging  behavior  of  the  double  layer  present  on  the  interface 
between  electrodes  and  electrolytes.  However,  such  models  do 
not  include  the  gas  and  temperature  dynamics  occurring  through 
the  flow  paths  and  in  the  cell  when  the  current  is  applied.  In  addi¬ 
tion,  the  partial  pressures  of  reactants  drop  along  the  pores  in  the 
gas  diffusion  layer  (GDL)  that  affect  the  net  pressure  exerted  on 
the  chemical  reaction  rate  and  increase  the  over-potential,  are 
not  considered.  The  temperature  rise  that  eases  water  removal, 
increases  the  chemical  reaction  and  subsequently  affects  the  out¬ 
put  voltage  of  a  cell  has  not  been  taken  into  account  in  this 
model. 

On  the  other  hand,  the  CFD-based  models  have  been  widely 
used  to  analyze  transport  mechanisms  of  the  mass  and  charge 
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Fig.  1.  A  schematic  diagram  of  a  PEM  fuel  cell  system. 


and  their  spatial  distributions  in  a  single  cell,  but  are  limited  in 
their  ability  to  represent  the  dynamic  characteristics  of  a  stack 
for  control  purposes.  In  addition,  exponential  growth  of  compu¬ 
tational  time  required  for  unsteady  analyses  impedes  application 
of  the  model  for  design  of  controls. 

Our  model  was  devised  on  the  basis  of  empirical  equations 
and  takes  three  additional  major  factors  into  account:  membrane 
water  balance,  diffusion  in  the  GDL  and  temperature  distribu¬ 
tion. 

A  PEM  fuel  cell  is  constructed  by  connecting  individual  mod¬ 
els  for  layers.  The  I-V  characteristic  is  the  difference  between 
the  open  circuit  voltage  and  the  over-potentials,  which  include 
the  ohmic  losses  in  the  membrane,  the  activation  over-potential 
in  the  catalyst  layer  on  the  cathode  side  and  the  concentra¬ 
tion  over-potential.  The  relationship  for  a  single  cell  may  be 
written  as  a  function  of  the  physical  parameters  like  the  reac¬ 
tant  partial  pressure,  temperature,  current  and  membrane  water 
content  [6].  The  output  characteristic  of  a  stack  is  assumed 
to  be  a  multiple  of  a  single  cell  characteristic  and  given  as 
follows. 

Tcell  =  E(p,  T )  —  Pact  (Pi  T,  i)  —  F  ohmic  (A  ^membr>  T) 

-Fconc (P,  T,  i)  (1) 

Vst  =  ftPcell 

The  dynamics  of  a  fuel  cell  system  involve  mass  flows  of  air 
and  water.  The  air  supplied  flows  through  the  gas  flow  channel 
and  the  GDL  before  reaching  the  catalyst  layer,  and  at  the  same 
time,  takes  up  water  from  the  humidifier.  Water  generated  in  the 
catalyst  layer  diffuses  through  the  membrane  where  protons  take 
water  from  the  anode  to  the  cathode  side.  The  heat  generated  by 


the  chemical  reaction  and  charge  transport  raises  the  tempera¬ 
ture  in  the  cell.  All  of  these  changes  affect  the  dynamic  behavior 
of  the  cell.  Fuel  cell  dynamics  can  be  improved  by  manipulat¬ 
ing  three  factors:  (1)  water  dynamics  in  the  membrane,  (2)  the 
partial  pressure  of  oxygen  drop  in  the  GDL  and  (3)  temperature 
variation. 

Membrane  water  content  determines  proton  conductivity. 
The  dynamics  of  water  content  are  described  by  the  electro- 
osmotic  driving  force  due  to  the  different  electrochemical 
potentials  at  the  anode  and  cathode,  and  the  diffusion  caused 
by  the  water  concentration  gradient  at  the  two  boundaries  [10]. 
Considering  the  water  mass  that  flows  at  the  boundaries  of 
the  membrane  layer,  the  dynamics  of  the  water  concentra¬ 
tion  in  the  membrane  can  be  improved  as  follows  [11],  where 
C  is  the  mass  concentration  (kgm-3),  M  is  the  mole  mass 
(kg  mol-1),  b  is  the  membrane  extension  coefficient  [10],  p  is 
the  membrane  dry  density  (kg  m-3)  and  Aceu  is  the  fuel  cell  area 

(m2)- 

^  _ _ Ch2o  ,mass  /Mh2Q _ 

Airy, membr/ -Vnembr  ^Cj^O.mass  /Mh20 

_ d  (  Ch20,  mass  cell  Giembr)  (2) 

PI  water,  membr  —  ] 

d  t 

=  Ikele,  membr,  an  fkele,  membr,  caT"  W cliff,  membr,  an  T"  Wdiff,  membr,  ca 

The  electro-osmotic  driving  force  created  by  the  different 
electrochemical  potential  at  the  anode  and  cathode  determines 
the  water  mass  flows  of  Wele, membr, an  and  Wele, membr, ca  at  the 
boundaries  of  the  membrane  layer.  In  addition,  the  diffusion 
caused  by  the  water  concentration  gradient  at  the  two  bound¬ 
aries  makes  up  the  mass  flows  of  Wdiff, membr, an  and  Wdiff, membr, ca- 
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Those  relationships  are  described  by  Eqs.  (3)— (5)  proposed  by 
Springer  [10]. 


n&  =  0.0029+embr  +  0.05/,memhr  -  3.4  x  1(T19 

i 

H/ele,membr,z‘  =  ^water^ celled,?  ~ 

F 


Wdiff,membr,i  —  ^water^  cell  D water 


(Pi  E’mid) 

finembr 


(3) 

(4) 

(5) 


where  the  diffusion  coefficient  Dwater  and  the  water  concen¬ 
tration  Q  are  calculated  from  the  empirical  equations  [10] 
expressed  as  a  function  of  membrane  water  content  Amembr. 


D 


water  : 


f^(^membr)  exP  (  2416  ( 


1 


303  Tmembr 


(6) 


fX^-membr)  —  < 


kt6 

10  ^(l  +  2(Amembr  —  3)) 
10-6(3  -  1.67(A.membr  -  3)) 
1.25  x  IQ-6 


2  >  A,membr 

3  >  7membr  ^  2 
4.5  >  A,membr  >  3 
7-membr  4.5 


(7) 


The  boundary  water  content  A*  is  a  function  of  water  activity 
ai,  which  is  calculated  from  the  water  vapor  partial  pressure 


'  0.043+17.81^—  39.85a2+36tf3 

1  >  0/  >  0 

ki  =  < 

14+1.4(0i  -  1) 

3  >  0/  >  1 

(8) 

16.8 

V 

at  >  3 

J 

di  —  J- 

> 

(9) 

The  reactant  entering  the  cell  diffuses  through  the  GDL 
before  reaching  the  catalyst  layer,  and  significantly  affects 
the  overall  dynamics  of  the  reactants.  This  diffusion  effect  is 
described  by  using  the  mass  continuity  and  the  Stefan-Maxwell 
Eq.  (10) 


£<y  dpi  8Ni 

+  — -  =  0 
RT  dt  3y 

£g  dpi  RT 

PcaDik 


( PiNk  -  pkNt ) 


(10) 


Hence,  i,ke(  1,3)  represents  the  sum  of  the  species  partial 
pressures,  where  p\  is  the  oxygen  partial  pressure,  P2=Psat(T) 
and  p^  are  the  water  vapor  and  the  nitrogen  partial  pressures, 
respectively,  and  the  diffusion  coefficients  of  Pc^ik  include  the 
cathode  pressure  of  /?ca.  The  parameters  sg  and  t  are  the  porosity 
and  the  tortuosity  of  the  GDL. 

If  a  cell  assembled  with  cubical  layers  is  isotropic  and  con¬ 
stant,  then  the  energy  conservation  equation  can  be  applied. 
Accordingly,  the  total  energy  change  in  a  controlled  volume  is 
equal  to  the  sum  of  the  energy  exchanged  at  the  boundaries  of  a 
control  volume  and  internal  energy  resources.  In  fact,  the  energy 
exchanged  at  the  boundaries  occurs  in  two  ways:  (a)  conduction 
across  the  cell  and  (b)  convection  occurring  between  bipolar 
plates  with  the  coolant,  reactants  and  water.  Thus,  the  thermal- 
dynamic  behavior  can  be  described  using  the  following  energy 


conservation  Eq.  (11)  [11]: 

^  drcv 

/  P Pi  C^mass^celflcv " 


dt 


=  WinCPj(Tin  -  Tcv)  +  Q 


mass  flow  in 

+  Gcond  +  Gsou 

conduction  heat  transfer  sources 


conv 

convection  heat  transfer 


(ID 


The  internal  energy  source  is  composed  of  the  entropy  loss 
and  the  chemical  energy  required  for  protons  to  overcome  the 
barrier  of  the  over-potentials  in  both  catalyst  layers  (12).  In  addi¬ 
tion,  other  source  terms  are  ohmic  losses  caused  by  transport  of 
electrons  and  protons  in  the  cell  [13]: 


Gsou  —  7 


f  TAs 
\  4f 


+  f’act  +  f^membr 


(12) 


where  As  is  equal  to  —326.36  (J  mol  1  K  ])  [14],  uact  is  given 
in  the  reference  [15]  and  Rmembr  is  the  membrane  resistance. 


2.2.  Air  supply  system 

The  air  supply  system  consists  of  four  subsystems;  an  air 
supplier,  a  humidifier  and  an  inlet  and  outlet  manifolds  with 
a  regulator  adjusting  the  pressure  at  the  stack.  It  functions  to 
continuously  replenish  the  air  to  the  fuel  cell  stack  as  the  load 
varies.  As  an  air  supplier,  a  compressor  or  a  blower  are  widely 
employed.  One  of  the  advantages  for  the  use  of  a  blower  is  the 
less  consumption  of  the  parasitic  power  than  that  of  the  com¬ 
pressor  and  consequently  overall  efficiency  of  the  system  can  be 
augmented  [16].  The  blower  is  connected  between  a  humidifier 
and  inlets  of  flow  channels  via  pipes.  The  humidifier  described 
here  is  simplified  as  the  ideal,  with  no  associated  dynamics 
or  energy  losses.  The  blower  is  usually  driven  by  an  electric 
motor.  The  dynamic  characteristics  of  the  blower  system  can  be 
described  by  the  sum  of  all  moments  of  inertia  of  the  motor  and 
the  impeller,  and  the  torque  produced  by  the  motor.  Hence,  the 
torque  produced  by  the  motor,  Tbi,m  (7),  is  a  function  of  the  stator 
resistance,  Rs,b\,m  (Ohm),  flux  linkage,  0bl,m  (V  s  rad-1)  and  the 
number  of  the  poles,  ftbi,m,pi>  stator  voltage,  Vbi,m  (V) 

[17]. 


da>b\  1  ( 

~dt~  _  4b\  Vbl,m  ~ 

3  7?bl,m, 
Jb\,m  —  0b\,m~  ~ 


fkbl  APbl  ^?bl,m 
0b]  Pamb^bl 


pi  ^bl,m,pl 
7?s,bl,m 


/  Nbl,m,pl  \ 

V  2  J 


^bl,m&>bl 

(13) 


where  co  is  the  angular  velocity  (rads-1),  J  is  the  rotational 
inertia  (kg  m2),  r]  is  the  efficiency,  p  is  the  pressure  (Pa)  and  p  is 
the  air  density  (kg  m-3).  The  air  blower  flow  rate  is  a  function  of 
the  angular  velocity  and  pressure,  and  the  efficiency  is  a  function 
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of  the  flow  rate  and  the  angular  velocity  [17]. 


WU=< 


a>bi(— 20.581  x  ( p *)2  -  1.4415  x  10-3  x  p* 

+4.1333  x  10-5),  p*  <  9  x  10-4  Pas2  rad-2 

otherwise,  ct>bi( —  1 .7973  x  p*  +  1.6409  x  10-3) 

(14) 


m  =  -2.8831  X  1013  X 


3 

+9.5115  x  108  x 


+  1.3087  x  104  x 


+  0.17945 


where  p*  is  (pca  -  pamb/o->h  )- 

The  blower  parameters  are  extracted  by  characteristic  data 
and  specifications  as  provided  by  PADT  (Phoenix  Analysis  & 
Design  Technologies)  [18],  which  includes  both  the  flow  param¬ 
eter  and  overall  efficiency  versus  the  head  parameter. 

Dynamic  characteristics  of  the  inlet  and  outlet  manifold  pres¬ 
sures  are  described  by  using  the  mass  conservation  equation. 


Hence,  Fr  denotes  the  frontal  area  (m2)  of  the  radiator,  and 
^rad, c, in  denotes  the  radiator  inlet  coolant  temperature  (K).  The 
electric  power  for  the  fan  can  be  calculated  according  to  a  ther¬ 
mal  dynamic  relationship  between  the  pressure  drop  and  the  air 
flow  rate  [13], 

Pfan  =  - ' - (WmCKTamh\P(,k-[lk)  -  1])  (20) 

^elec^fan 

where  Pfan  denotes  the  electric  power  (W)  of  the  fan. 

The  reservoir  stores  the  coolants  after  a  heat  exchange 
at  the  radiator  by  convection.  The  variation  of  the  heat  in 
the  reservoir  is  the  sum  of  the  heat  that  the  coolant  car¬ 
ries  and  the  heat  exchanged  with  the  ambient.  Therefore,  the 
reservoir  outlet  coolant  temperature  at  the  end  of  the  given 
time  interval,  rres,c,out  (K),  can  be  expressed  by  the  equation 
[13] 

^  (rr 

Tres,c,out  =  Tj-es,p  —  - —  (WcC„(Tres,p  -  Ties,  c,  in) 

^res  Cpres 

+  ^res2lres(Tres,c,in  -  Tamb))  (21) 


=  ^(WbiTbi  -  Wim,out7im)  (16) 

At  Vim 

2.3.  Thermal  circuit 

A  thermal  circuit  should  dissipate  the  heat  produced  by  the 
stack.  The  circuit  has  3  components  including  a  three-way  valve 
to  allow  the  coolant  to  bypass,  a  radiator  to  exchange  heat  with 
the  ambient  media,  a  fan  to  increase  the  effectiveness  of  heat 
convection  and  a  reservoir  to  store  the  coolants.  A  coolant  pump 
supplies  the  coolant  for  the  heat  source. 

The  opening  width  of  the  bypass  valve  is  assumed  to  be  linear 
with  a  factor  k.  Then,  the  coolant  temperature  at  the  reservoir 
inlet,  rres, c, out,  can  be  expressed  as  a  function  of  the  stack  outlet, 
Tit, c, out,  and  the  radiator  outlet,  Tra( j,c,0ut  and  k. 

WcCpc  Ties, c, in  =  (1 — k)WcCpc  Tsi  c  oui~\-kWcCpc  Tlad,c,out  (IT) 

The  behavior  of  the  radiator  follows  thermodynamic  princi¬ 
ples  Kroger  [19]  proposed  an  empirical  equation  for  the  heat 
transfer  coefficient  of  the  radiator,  /zrac i  (kWm-2  K-1)  and  pres¬ 
sure  drop,  pv  (kPa)  of  the  radiator  as  a  function  of  the  air  flow 
rate,  Wair  (kgs-1)  [13]. 

hrad  =  —1.4495 W2r  +  5.9045  l+ajr  -  0.1157 

Pr  =  (326. 12 Vkair  -  75.396)  +  101.325  ' 

If  the  heat  of  the  coolant  is  transferred  completely  to  the  radi¬ 
ator  without  any  loss,  the  heat  capacity  of  the  coolant  is  identical 
to  that  of  the  radiator.  Thus,  the  radiator  outlet  coolant  temper¬ 
ature  can  be  expressed  as  a  function  of  the  radiator  geometry 
and  the  heat  convection  caused  by  the  temperature  difference 
between  the  ambient  and  the  radiator  outgoing  air  temperature 

[13]. 

rr  rj,  a  s'  ( Fr(Tlad,c,in  —  Tamb)/lrad  \  /inx 

Tlad,c,out  =  Tlad,c,in  -  0.46  - —— -  (19) 

V  wccPc  j 


where  Tles,p  is  the  temperature  (K)  of  the  reservoir  at  the  previ¬ 
ous  of  time  point,  At  is  the  time  interval  (s),  mres  is  equivalent 
to  the  coolant  mass  in  the  reservoir  (kg),  rres,c,in  is  the  reser¬ 
voir  inlet  coolant  temperature  (K)  and  hVQS  is  the  convective 
heat  transfer  coefficient  of  reservoir  to  the  ambient  temperature 
(W  m-2  K-1),  respectively. 

Assuming  that  all  of  the  heat  generated  in  the  stack 
is  completely  transferred  to  the  coolant;  the  coolant  flow 
rate  is  expressed  by  its  relationship  with  the  heat  source 

[13]. 


Wc  = 


Gsou 

Cpc  AT 


(22) 


At  equilibrium,  the  excessive  heat  rejected  by  the  coolants 
is  identical  to  the  sum  of  the  heat  the  reservoir  stores  and 
that  which  the  radiator  exchanges  with  the  ambient  environ¬ 
ment.  If  a  maximum  coolant  flow  rate  is  assumed  then  the 
temperature  drop  along  the  coolant  flow  channel  can  be  cal¬ 
culated  by  considering  that  the  maximum  heat  produced  in  the 
stack  should  be  the  same  as  the  amount  of  heat  carried  by  the 
coolant  flow  rate.  The  temperature  at  the  outlet  of  the  coolant  is 
assumed  to  be  345.15  K  because  the  temperature  gradient  from 
the  catalyst  to  the  coolants  channel  is  8  K  at  a  maximum  load 
current  if  the  catalyst  temperature  is  353.15  K.  Thus,  the  inlet 
temperature  is  obtained  at  a  temperature  drop  of  12  K  and  a 
maximum  flow  rate  of  3 kgs-1.  The  heat  transfer  coefficient 
of  the  radiator  and  the  volume  of  the  reservoir  were  chosen 
on  the  basis  of  the  maximum  heat  capacity  stored  and  dissi¬ 
pated. 


3.  Control  strategies 

Fig.  2  depicts  a  block  diagram  of  the  air  supply  system  and 
thermal  circuit  with  either  PI  controllers  or  state  feedback  con¬ 
troller. 
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Fig.  2.  Block  diagram  for  an  air  and  coolant  control  with  (a)  classic  PI  controllers  and  (b)  state  feedback  controller  with  integral  controller. 


The  feedback  loop  of  the  air  supply  system  should  main¬ 
tain  the  optimal  oxygen  ratio  and  prevent  oxygen  starvation  that 
might  occur  during  abrupt  changes  in  the  load  current.  Hence, 
the  oxygen  excess  ratio  is  defined  as  the  rate  of  oxygen  sup¬ 
plied  to  that  consumed.  While  the  amount  of  oxygen  consumed 
depends  on  the  stack  current,  the  amount  of  oxygen  supplied  to 
a  fuel  cell  is  directly  related  to  the  blow  motor  voltage.  Thus, 


the  air  supply  controller,  a  static  feed-forward  controller  (sFF) 
[20],  uses  a  polynomial  that  interpolates  map  data  and  includes 
an  optimal  relationship  between  the  stack  current  required  and 
the  motor  voltage  of  the  blower  in  order  to  maintain  the  oxygen 
ratio  at  2.  In  this  case,  the  air  flow  of  the  blower  can  simply  be 
controlled  by  the  blower  motor  voltage.  The  design  of  this  type 
of  controller  has  extensively  reviewed  by  other  authors  [21-23]. 
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As  shown  in  Fig.  6,  the  rejection  behavior  of  the  oxygen  excess 
ratio  at  an  occurrence  of  a  temperature  rise  can  be  optimized  by 
the  controllers. 

The  stack  is  regarded  as  a  single  thermal  mass  with  a  heat 
capacity.  Under  the  assumption  that  the  heat  exchange  by  radi¬ 
ation  and  convection  to  the  ambient  is  negligible,  and  the  stack 
temperature  is  equal  to  the  average  of  the  stack  outlet  coolant 
temperatures  on  the  anode  and  cathode  side,  the  variation  of  the 
temperature  in  the  stack  is  equal  to  the  sum  of  the  heat  source 
terms  in  the  stack,  the  heat  exchanged  with  the  coolants. 

mstCpst  ^  =  <2 sou  H-  WcC pc(Tc,in  —  Ts t)  (23) 

where  mstCPst  is  the  heat  capacity  of  the  stack  (J  K-1),  Wc  is  the 

coolant  flow  rate  (kgs-1),  the  control  variable  and  2S0U  is  the 
internal  energy  source  (J  s-1),  a  function  of  the  load  current. 

Due  to  the  nonlinearity  of  the  equations  in  the  lumped 
thermal  stack  (Eq.  (23))  and  the  reservoir  model  (Eq.  (21)), 
Talyor’s  expansion  is  used  for  linearization  at  an  operating  point, 
where  the  reservoir  temperature  and  coolant  flow  rate  are  set  at 
337.15 K  and  0.96 kgs-1,  respectively.  The  stack  current  and 
voltage  are  140  A  and  198  V,  respectively.  The  state  equations 
and  variables  are  defined  as  follows: 

8x  ==  A8x  “I-  Buu  “I-  B™w 

(24) 

8y  =  C8x 

x=  [  Tst  Tre  s  ]  (states) 
u  —  [Wc  k  ]  (controlled  input) 
w  =  7st  (disturbance) 
y=  [Tst  Ties  ]  (output) 


Gci(s)  =  Kp,  c  +  —  (26) 

Eq.  (21)  describes  the  relationship  between  the  tem¬ 
perature  in  the  reservoir  and  the  factor  for  the  bypass 
valve  opening,  which  the  transfer  function  is  given  in 
Eq.  (27).  Likewise,  the  gains  for  the  PI  controller  are 
selected  so  that  the  bandwidth  of  the  closed-loop  is  5  times 
higher  than  the  time  constant  of  the  coolant  flow  feed¬ 
back  outer  loop.  In  addition,  the  damping  ratio  is  set  to 
0.707.  The  resulting  gains  are  Kv b  =  0.1902  and  Ki^  =0.0546 
(s-1). 

Tresis)  =  (W^Cpc  +  ^res^resXTamb  T^)  (27) 

k(s)  ( mIesCpresS  +  W^Cpc) 

GC2(S)  =  Kpt>  +  ^  (28) 

3.2.  Design  of  state  feedback  controls  with  integral  controls 

The  parasitic  power  dissipated  in  the  coolant  pump  is  not 
considered  as  a  control  object  by  designing  the  classic  PI  con¬ 
troller,  even  though  it  can  reject  sufficient  amounts  of  heat 
and  effectively  suppresses  temperature  surges  in  the  layers  of 
the  cell.  One  alternative  is  the  use  of  a  state  feedback  con¬ 
trol,  where  the  parasitic  power  dissipated  in  the  coolant  pump 
can  be  considered  to  be  one  of  the  control  objectives.  On 
the  other  hand,  the  parasitic  power  of  the  coolant  pump  is 
directly  proportional  to  the  coolant  flow  rate.  Thus,  the  coolant 
flow  rate  is  included  as  a  parameter  in  the  cost  function  given 
below.  Optimization  of  the  gains  is  achieved  by  the  LQR  (linear 
quadratic  regulator)  method,  which  sums  the  square  of  the  errors 
[24]. 


where  the  matrices  of  the  linearized  system,  A,  BU,BW  and  C  are 
listed  in  the  Appendix  A. 


J  = 


(8xT  Qx8x  +  8uT  R8u)dt 


(29) 


3.1.  Design  of  classic  PI  controls 

The  state  equations  derived  above  present  a  multi-input- 
multi-output  structure,  where  two  controlled  input  variables, 
coolant  flow  rate  and  bypass  valve  opening  factor,  are  depen¬ 
dent  on  each  other.  This  dependence  can  be  minimized  if  the 
time  constants  of  two  feedback  loops  are  set  in  a  different  order. 
In  that  case,  the  temperature  in  the  stack  can  be  controlled  by  the 
coolant  flow  rate  independent  of  the  temperature  of  the  coolant 
being  controlled  by  the  opening  factor  k.  Eq.  (23)  includes  a 
relationship  between  the  stack  temperature  and  the  coolant  flow 
rate,  whose  transfer  function  is  given  in  Eq.  (25).  The  system 
follows  the  first  order  of  an  ordinary  differential  equation  and 
thus  a  classic  PI  controller  is  employed.  The  two  gains  of  the  PI 
controller  are  selected  by  the  bandwidth  of  the  closed-loop  that 
is  3  times  higher  than  the  time  constant  of  the  heat  source  term 
response  and  a  damping  ratio  of  0.707.  The  resulting  gains  are 
Kpx  =  0.25  and  Klc  =  0.017  (s_1). 

W  =  Cp c  •  (T°s  -  t°) 

WcCs)  (mstCpst  ■  s  +  W°  ■  Cpc) 


Qx  represents  the  weighting  matrix,  which  amplifies  the  errors 
of  the  control  objects,  while  the  other  weighting  matrix  R  is  used 
to  suppress  the  effect  of  manipulating  the  variables. 

The  state  equation  of  the  control  plant  presents  a  2-by-2 
matrix,  where  the  variables  are  coupled  with  each  other.  A 
decoupling  of  the  two  loops  has  been  accomplished  by  assigning 
different  time  constants  to  the  two  closed  loops.  In  fact,  the  valve 
opening  factor  does  not  directly  affect  the  dynamics  of  the  stack 
temperature,  while  the  reservoir  temperature  is  strongly  influ¬ 
enced  by  the  valve  opening  factor  rather  than  the  coolant  flow 
rate.  Hence,  the  time  constant  of  the  transfer  function  between 
the  stack  temperature  and  coolant  flow  rate  is  5  times  faster 
than  that  between  the  stack  temperature  and  the  valve  opening 
factor. 

On  the  other  hand,  integrators  are  required  to  suppress  any 
steady  state  errors.  Thus,  the  errors  of  both  closed  loops  are 
defined  as  a  new  state  variable  that  is  considered  in  the  cost 
function; 


(30) 
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poo 

J=  (8xT  Qx8x  +  qT  Qiq -\- 8uT R8u)dt  (31) 

Jo 

where  Q\  is  the  weighting  matrix  for  integrator. 

The  rules  for  the  optimal  control  inputs  are  obtained  [25]; 

8u  =  -K[8x  q]T  =  -Kp-8x-  Ki-q  (32) 

where  the  controller  gain  is  K=R~l Bu/TP.  P  is  the  solution  of 
the  Algebraic  Riccati  equation  that  is  given  as  follows: 

PA' +  A' P  +  Q  -  PB^R-1  B'l P  =  0  (33) 

where  A\  B'u,  Q  =  dia g(Qx,  Q\)  and  R  is  listed  in  Appendix  A. 

When  the  weighting  matrix  R  is  larger  than  the  weighting 
matrix  Q,  the  role  of  the  coolant  flow  rate  in  the  cost  func¬ 
tion  increases,  and  subsequently,  the  gains  of  the  controller  are 
chosen,  which  minimizes  the  parasitic  power.  After  several  itera¬ 
tions  with  different  weighting  factors,  the  optimal  control  matrix 
Kv  and  K\  is  given  by 


—  1.2015 

0.01 

,  Ki= 

-0.1  0.0027" 

-0.0521 

-3.1479 

-0.0027  -0.1 

( 

3.3.  Disturbance  compensation 

The  heat  produced  in  the  stack  tends  to  follow  the  current 
drawn  from  the  stack.  The  current-dependent  heat  is  regarded  as 
a  disturbance  in  the  control  loop.  This  cannot  be  fully  rejected 
by  the  typical  coolant  control,  which  measures  the  temperature 
at  the  outlet  of  the  coolants.  As  a  result,  the  heat  dissipated  is 
less  than  the  heat  produced.  A  countermeasure  is  to  estimate  the 
temperature  rise  in  a  layer  that  is  directly  related  to  the  magnitude 
of  the  current  load  and  feed-forward  it  to  the  temperature  control 
loop  shown  in  Fig.  2.  The  relationship  between  the  current  and 
the  stack  temperature  is  derived  by  using  Eqs.  (12)  and  (23) 
yielding  the  following  transfer  function  (Eq.  (35)): 

tJa)  _  (~tas/4f  +  v^) 

I(s)  (mstCpst  ■  s  +  CpcW®) 

TJs)  =  flmembr  ^  ’ 

I2(s)  (mstCpst  ■  s  +  CpcW?) 

However,  all  previously  published  air  supply  control  designs 
assumed  a  constant  stack  working  temperature  [20-23].  In  fact, 
distribution  of  the  temperatures  in  the  individual  layers  through 
the  plane  varies  because  of  the  various  sources  of  energy  occur¬ 
ring  in  the  chemical  reactions,  and  Joule’s  losses  associated  with 
charge  transport.  In  order  to  dissipate  the  heat  in  the  stack,  the 
temperature  of  the  coolant  control  loop  is  set  lower  than  that  in 
the  stack,  thereby  changing  the  temperature  in  the  gas  flow  chan¬ 
nels.  When  the  temperature  in  the  channel  is  lower,  the  pressure 
drops  according  to  the  ideal  gas  law  for  the  given  volume,  and 
then  a  pressure  difference  to  the  inlet  manifold  becomes  larger. 
On  the  other  hand,  the  mass  flow  rate  at  the  inlet  of  the  stack 
increases  according  to  the  nozzle  equation  [6]  and  the  oxygen 
excess  ratio  increases.  The  surplus  air  is  reduced  by  an  additional 


element  in  the  controls  that  depends  upon  the  coolant  temper¬ 
ature  in  addition  to  the  current  that  determines  the  amount  of 
oxygen  consumed.  However,  the  relationship  between  the  load 
current  and  temperatures  given  as  inputs  and  the  blower  voltage 
as  an  output  is  nonlinear  at  an  optimum  oxygen  excess  ratio. 
Therefore,  the  set  of  data  for  the  relation  is  obtained  from  multi¬ 
ple  runs  of  the  entire  model  at  different  currents  and  temperatures 
and  used  for  a  compensation  of  the  effects. 


4.  Simulation  and  analyses 

Simulations  are  performed  to  analyze  the  dynamic  behav¬ 
ior  of  the  stack  along  with  the  air  supply,  thermal  system 
and  the  associated  control  strategies.  Dynamics  of  water  con¬ 
tent  in  the  membranes,  temperature  variations,  oxygen  excess 
ratio  and  system  responses  at  load  currents  are  discussed 
in  the  following  sections.  The  empirical  parameters  and  ref¬ 
erence  data  for  the  models  chosen  are  given  in  Table  1 
([6,10,12,13,26]).  All  models  were  coded  by  blocks  given  in 
MATLAB/Simulink. 
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Table  1 

Simulation  parameters 


Fuel  cell 

n 

381 

Me 

0.028  m2 

Proton  conducting  model  [12] 

*1! 

0.5139 

b\2 

0.326 

b2 

350 

rid 

/(C water) 

Av 

/(T  C  water) 

Gas  transport  model  [10] 

Deff 

f(P,  T)m2  s-1 

P  sat 

AT) 

Electrochemical  reaction  model  [6] 

Po 

1.0  bar 

Tref 

353.15  K 

Eref 

1.229  V 

^catl,eff/ACell 

AT  A  Po2 ) 

Thermal  model  [13] 

Hgas 

Ap,T) 

C  Pgas 

AP,T) 

Pgas 

f(P  ’  T) 

Fr 

0.6 

Wires 

5 

brad 

5.897 

Thickness  (m)  Density  (kg  m  3)  Heat  conductivity  (W  m  lK  x)  Specific  heat  (J  kg  1  K  x) 


Geometrical  data  for  layers  [26] 


Coolant  channel 

0.001 

1400 

52 

935 

Plate 

0.001 

1400 

52 

935 

Gas  channel 

0.001 

1400 

52 

935 

GDL 

0.0004 

2000 

65 

840 

Catalyst  layer 

0.000065 

387 

0.2 

770 

Membrane  layer 

0.000183 

1967 

0.21 

1100 

4.1.  Water  content  in  the  membrane 

Fig.  3  shows  a  comparison  of  the  membrane  water  content 
between  the  empirical  and  the  proposed  model  at  a  step  load 


current.  The  membrane  water  content  is  a  function  of  the  water 
activity  (Eq.  (9))  that  depends  upon  the  saturated  vapor  pressure 
and  vapor  pressures  of  the  cathode  and  anode  sides.  Hence,  the 
saturated  vapor  pressure  is  a  function  of  the  temperature.  Since 
the  empirical  model  assumes  a  constant  temperature  of  353. 15  K 


Fig.  4.  Temperatures  of  the  catalyst  layer  and  coolant  channel  by  the  coolant  Fig.  5.  Temperature  variation  in  a  cell  depending  on  currents  with  and  without 
flow  controls.  the  feed-forward  of  the  disturbance. 
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Fig.  6.  Comparison  of  oxygen  excess  ratio  at  a  constant  and  varying  temperature. 
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Fig.  7.  Comparison  of  the  oxygen  excess  ratio  before  and  after  temperature 
compensation. 
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Fig.  8.  Comparison  of  (a)  coolant  flow  rate  and  (b)  stack  inlet  coolant  temperature  with  a  given  current  step  between  PI  controls  and  state  feedback  controls. 
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in  the  membrane,  the  vapor  pressure  likely  correlates  with  the 
change  of  the  load  current,  but  no  water  transfer  dynamics  are 
involved. 

Conversely,  the  variation  of  temperature  and  the  change  of 
water  balance  in  the  membrane  dynamically  influence  the  water 
content  in  the  membrane.  When  the  temperature  in  the  coolant 
channels  is  maintained  at  349.15  K,  the  temperature  in  the  gas 
flow  channel  is  higher  than  349.15  K  and  the  saturated  vapor 
pressures  in  both  of  gas  channels  get  increased.  Consequently, 
the  water  activity  and  the  water  content  become  lower.  The 
change  of  water  balance  is  strongly  influenced  by  the  electro- 
osmotic  force  that  is  a  function  of  the  current  (Eq.  (4)).  At  a 
high  load  current,  the  elevated  temperature  of  the  stack  leads 
to  a  high  saturated  vapor  pressure  and  a  low  relative  humidity 
on  both  sides  of  the  cell  as  well  as  high  electro-osmotic  driving 
force.  As  a  result,  the  membrane  water  content  gets  lower. 

4.2.  Temperature  in  the  cell 

Fig.  4  shows  the  temperature  in  the  catalyst  and  coolant 
channel  with  and  without  the  feed-forward  of  the  disturbance. 


The  actual  temperature  of  the  stack  is  usually  measured  at  the 
stack  outlet  coolants  on  the  anode  and  cathode  sides  and  aver¬ 
aged.  A  reference  temperature  for  the  coolant  control  is  set  at 
349.15  L. 

When  a  multi-step  current  is  applied  to  the  stack,  the  tem¬ 
perature  in  the  stack  rapidly  rises,  particularly  in  the  catalyst 
on  the  cathode  side.  The  temperature  rise  is  3-7  K  higher  than 
the  average  temperature  in  the  stack,  where  the  coolant  temper¬ 
ature  is  fully  controlled  at  the  reference  temperature  349.15  K 
(see  the  dotted  line  in  Fig.  4).  It  should  be  noted  that  the  cat¬ 
alyst  and  membrane  layers  could  be  overheated  and  thereby 
damaged. 

The  difference  in  temperature  in  the  layers  can  be  reduced 
by  a  feed-forward  (FF)  of  the  disturbance  to  the  coolant  control 
loop  that  should  reject  excessive  heat  as  quickly  as  possible. 
The  transfer  function  of  the  disturbance  is  given  in  Eq.  (35). 
The  result  of  the  control  strategy  proposed  is  illustrated  in  the 
Fig.  4  as  a  straight  line,  where  the  temperature  of  the  catalyst 
layer  is  maintained  close  to  353.15  K.  The  coolant  temperature 
shows  the  variation  in  the  catalyst  temperature.  However,  an 
instant  rise  in  temperature  cannot  be  fully  suppressed  because 


Fig.  9.  (a)  FUDS  and  a  current  profile  with  a  base  load  of  120  A,  (b)  temperature  of  the  catalysts  and  coolants  without  the  FF,  (c)  temperature  of  the  catalysts  and 
coolants  with  the  FF,  and  (d)  oxygen  excess  ratio  after  a  temperature  compensation. 
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of  the  high  thermal  mass  and  large  heat  capacity  of  the  stack.  In 
addition,  there  is  a  steady  state  error  caused  by  the  temperature 
difference  between  the  coolant  channel  measured  and  the 
catalyst  layer.  Nevertheless,  the  cell  is  cooled  effectively  and 
the  exposure  to  heat  on  individual  layers  is  minimized. 

Fig.  5  shows  the  effects  of  the  coolant  controls  on  temper¬ 
ature  distributions  through  the  plane  of  a  cell.  As  the  amplitude 
of  the  current  changes  stepwise  from  0.5  A  to  0.55  A,  0.65  A 
and  0.7  A,  the  stack  temperature  becomes  higher.  When  the 
feed-forward  is  applied,  the  overall  stack  temperature  is  lowered 
and  the  catalyst  temperature  is  maintained  at  353.5  K,  which 
is  significantly  lower  than  before.  Likewise,  the  maximum 
difference  of  the  temperature  between  the  catalyst  on  the 
cathode  side  and  the  coolant  channel  becomes  4  K  lower  than 
before.  As  a  result,  the  cooling  of  the  stack  is  more  effective. 

4.3.  Oxygen  excess  ratio 

Fig.  6  illustrates  the  oxygen  excess  ratio  at  constant  and 
dynamically  varying  temperatures  with  a  coolant  flow  control. 
Because  of  the  change  of  pressure  in  the  gas  flow  chan¬ 
nel  that  is  caused  by  variations  in  the  stack  temperature,  the 
oxygen  excess  ratio  is  inversely  changed  as  to  the  current 
change. 

Fig.  7  shows  a  comparison  of  the  oxygen  excess  ratio  before 
and  after  compensation  of  the  temperature  influence  on  the  air 
control  loop.  The  compensation  enables  the  oxygen  excess  ratio 
to  be  maintained  at  level  2,  even  though  the  current  applied  to 
the  stack  varies  stepwise,  which  implies  that  the  parasitic  power 
at  the  blower  is  less. 

4.4.  Comparison  of  the  PI  and  state  feedback  controls 

Comparison  between  both  controls  shows  that  the  parasitic 
power  of  the  state  feedback  controls  a  multi-step  current  load 
is  5%  less  than  that  of  the  PI  controller.  However,  the  dynamic 
response  is  improved  by  the  state  feedback  controls.  In  Fig.  8, 
step  responses  of  two  controls  are  simulated  with  the  models 
discussed  earlier.  The  output  states  are  the  coolant  flow  rate  and 
the  stack  inlet  coolant  temperature.  The  rise  time  of  the  coolant 
flow  rate  by  the  state  feedback  controller,  6  s,  is  4  times  faster 
than  that  found  with  one  PI  controls.  Likewise,  the  rise  time  of 
the  stack  inlet  coolant  temperature  is  3  times  faster  than  that  of 
the  PI  controller. 

Parasitic  power  is  calculated  as  the  sum  of  the  electrical  power 
necessary  for  driving  the  blower  and  the  coolant  pump.  The  con¬ 
trol  strategies  with  the  proposed  state  feedback  control  require 
100  kWs  at  the  multiple  step  current,  while  that  for  the  PI  control 
loop  requires  106  kWs. 

The  response  of  the  state  feedback  control  was  compared 
with  a  conventional  one  by  using  a  current  profile  obtained  from 
a  vehicle  tested  at  the  FUDS.  Fig.  9  shows  the  simulation  results 
for  the  two  different  control  strategies  at  the  load  current.  The 
peak  temperature  in  the  catalyst  layer  is  6  K  higher  than  the  work¬ 
ing  stack  temperature  by  the  control  without  the  feed-forward 
(FF)  control,  even  though  the  coolants  are  controlled  around  the 
set  reference  temperature  shown  in  Fig.  9b.  Fig.  9c  shows  the 


temperature  of  the  catalyst  and  coolants  with  the  FF  of  the  dis¬ 
turbance.  The  peak  of  the  temperature  is  similar  to  that  of  the 
others  for  the  first  200  s,  but  is  suppressed  with  a  time  delay 
associated  with  the  closed-loop  dynamic  compare  to  Fig.  9b. 
The  excursion  duration  of  the  catalyst  temperature  is  decreased, 
and  the  heat  energy  imposed  on  the  thin  layers  is  reduced,  which 
significantly  reduces  the  heat  stress  on  the  layers.  Correspond¬ 
ingly,  the  oxygen  excess  ratio  shown  in  Fig.  9d  is  maintained  at 
the  optimum  value  compared  to  Fig.  6. 

5.  Conclusion 

This  article  addresses  the  design  of  temperature  control 
strategies  for  PEM  fuel  cells  and  their  effects  on  dynamics  and 
performance.  We  used  a  dynamic  stack  model  and  considered 
gas  diffusion  in  the  GDL,  dynamic  water  balance  in  the  mem¬ 
brane,  temperature  variation,  and  components  of  the  air  supply 
and  thermal  system.  There  are  several  major  outcomes: 

•  Dynamic  stack  behavior  was  improved  by  adding  dynamic 
water  balance  in  the  membrane  and  a  partial  pressure  drop 
of  reactants  in  the  GDL  and  temperature  distributions.  The 
results  show  that  distribution  of  the  temperature  through 
the  plane  is  asymmetric  and  the  temperature  rise  amounts 
to  3-7  K,  which  has  the  potential  to  damage  the  layers  at 
a  high  current  load.  Therefore,  proper  control  of  the  air 
and  temperature  might  be  required  to  ensure  durability  and 
efficiency; 

•  Most  control  strategies  focused  on  optimization  of  the  air 
supply  system,  where  the  working  temperature  in  the  fuel 
cell  stack  is  presumed  to  be  constant.  However,  we  find  that 
the  oxygen  excess  ratio  varies  inversely  with  temperature. 
Thus,  the  ideal  oxygen  excess  ratio  necessary  to  prevent  oxy¬ 
gen  starvation  cannot  be  maintained  at  the  optimum  value 
of  2; 

•  New  control  strategies  proposed  include  a  state  feedback  con¬ 
trol  with  a  feed-forward  of  the  disturbance  and  a  compensator 
for  minimization  of  the  temperature  effect  on  the  air  flow  rate. 
For  the  design  of  the  temperature  controller,  the  thermal  cir¬ 
cuit  is  approximated  with  a  second  order  system.  Classic  PI 
and  state  feedback  controls  are  used  to  compare  the  effective¬ 
ness  of  cooling.  The  results  show  that  the  temperature  rise  in 
the  catalyst  can  be  kept  within  an  allowable  value  and  dura¬ 
tion.  In  addition,  the  oxygen  excess  ratio  can  be  maintained 
at  an  optimal  value  by  minimizing  the  influence  of  temper¬ 
ature  variations  in  the  gas  flow  channel.  Consequently,  the 
power  consumption  of  the  blower  can  be  reduced  by  more 
than  15%  by  compensation,  and  5%  by  the  controlling  bypass 
valve  at  a  multi-step  load  profile.  Using  these  techniques, 
the  total  parasitic  power  was  reduced  by  approximately 
7%. 

Future  work  will  assess  other  factors:  (1)  design  of  an 
observer  for  the  temperature  in  the  layers  aiming  at  advanced 
controls  and  real-time  diagnosis  as  well  as  the  water  content  in 
the  membrane  and  (2)  optimization  of  the  air  and  temperature 
controls  with  a  real  behavior  of  a  humidifier. 
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